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The spin- wave theory and the coherent potential approximation are applied to a spin S Heisen- 
berg antiferromagnet with nonmagnetic impurities on square lattice. The impurity effects are 
taken into account by substituting 5(1 — x) for S and using the coherent potential approximation 
to the exchange interaction, where x is the impurity concentration. At T — for S — 1/2 the 
critical impurity concentration x c of the Neel state is 0.303 and the percolation threshold x p is 
0.500. The ground state in x c < x < x p is the disordered state with the spin gap. For S > 1 
the long range Neel order vanishes at x p — 0.500. These results explain qualitatively the experi- 
mental results of La 2 Cui_ a: Mg :E 04 (5 = 1/2) and K 2 Mm_ x Mg :I .F4 (5 = 5/2). The difference of 
x c between these materials is caused by the decrease in the magnitude of the effective spin with 
impurity doping. The spin gap is expected to be observed for L^Cui-rrMg^CU in x c < x < x p 
at low temperatures. 
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§1. Introduction 

Two dimensional magnet with nonmagnetic impuri- 
ties shows interesting properties in relation to the mag- 
netic order and the percolation. [1-12] The Neel temper- 
ature of an antiferromagnet is depressed with nonmag- 
netic impurity doping. The critical impurity concentra- 
tion for the two dimensional (2D) Ising antiferromagnet, 
e.g., K2Coi- 2; Mg ;:! .F4 agrees well with the percolation 
threshold 0.41. The critical concentration for the 2D 
Heisenberg antiferromagnet, e.g., K2Mni_ x Mg :E F4 also 
agrees with the percolation threshold 0.41. The agree- 
ment between the critical concentration and the perco- 
lation threshold can be understood by the fact that the 
long range order cannot exist if the percolation path docs 
not connect to infinity. In compared with above two ma- 
terials, by extrapolating the experimental data it is found 
that the critical concentration of L^Cui-zMg^O^ ex- 
pressed by the antiferromagnetic Heisenberg model with 
nonmagnetic impurities on square lattice, is not the per- 
colation threshold 0.41 but 0.2-0.25. [1] We consider this 
difference may be caused by the magnitude of the spin, 
5, of the host material, i.e., Cu 2+ (5 — 1/2) of La2Cu04 
and Mn 2+ (5 = 5/2) of K 2 MnF 4 . 

If the spin operator of the Heisenberg model on square 
lattice is approximated by the classical vector, we can 
obtain the spin dependence of the critical concentration 
of the Neel state. In this approximation, the Neel state 
energy is —ANJS 2 and the singlet dimer state energy is 
—NJS(S+1), where N is the lattice number and J is the 

* Present address : Institute for Solid State Physics, The Univer- 
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exchange integral. If the nonmagnetic impurity is doped 
and we adopt the coherent potential approximation, then 
J would be replaced by the coherent exchange integral 
J p . Comparing these energies, we find that the Neel state 
is stable when 5 > 1/3. However, spin 5 is allowed only 
to be an integer or a half odd number. This l/3(= 5 crl ) 
may be regarded as the effective spin S eS averaged over 
all impurity configurations. As the spin 5 is reduced 
to 5 crl by the impurity doping, the critical impurity con- 
centration x c satisfies the relation, 5(1 — x c ) — 5 crl = |. 
At 5 — 1/2 this critical concentration x c is about 0.3 
from ^(1 — %c) — Thus if the nonmagnetic impurity 
doping reduces the effective spin, the singlet state can be 
stable when S oS < 1/3 and x > 0.3. On the other hand, 
at 5 = 5/2 the critical concentration x c is about 0.87 
and is larger than the percolation threshold 0.41. Thus 
the antiferromagnetic order for 5 = 5/2 persists up to 
the percolation threshold. 

The critical concentration for the antiferromagnetic 
Heisenberg model with 5 = 1/2 on square lattice 
has been investigated by the present authors using the 
Kondo-Yamaji's Green's function method with nonmag- 
netic impurities. [13] The obtained results show that the 
critical concentration is not equal to the percolation 
threshold and the disordered state which is character- 
ized by the finite spin gap and the exponentially decay- 
ing correlation function appears in the range of x c < x < 
x p . However, the generalization of the Kondo-Yamaji's 
method to the general spin 5 is difficult, because we 
cannot self-consistently determine a new decoupling pa- 
rameter. On the other hand, the spin-wave theory can 
be applied to any spin systems. Therefore, in this paper 
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we study by using the spin-wave theory and the modified 
spin- wave theory. [14] 

The present paper is organized as follows: in § 2, we 
formulate the spin-wave theory for a 2D Heisenberg an- 
tiferromagnet with nonmagnetic impurities and obtain 
the critical concentration and the percolation threshold 
for any spin 5 at T = 0; in § 3, the modified spin-wave 
theory is applied in the range between the critical con- 
centration and the percolation threshold. Section 4 is 
devoted to summary and discussion. 

§2. Spin-wave theory of a 2D Heisenberg anti- 
ferromagnet with nonmagnetic impurities 

We consider the antiferromagnetic (J > 0) Heisenberg 
model with nonmagnetic impurities on square lattice: 



H c = 2J y OiOjSi ■ Sj , 



<i,3> 



(1) 



where X)<ij> denotes the summation over all nearest 
neighbor sites. The magnetic occupation operator Oi are 
defined by 

1 if the site i is occupied by a spin , 
otherwise . 

(2) 

The subscript c stands for a given spin configuration. 
Assuming the lattice to be bipartite, we use i, 1,1' 6 A 
sublattice and j, m,m' 6 B sublattice. We introduce the 
antiferromagnetic Holstcin-Primakoff transformation for 
the impure case: 



aiS, =aia}\2S 



aia\ai 



er;5 ; + = \J2S' - (Jiajai am 
aiSf = 5' - majai 



<J rnS m 



2S" 



a m S+ = a m bL\ 2S 



&m &m &m. 



h 1 h 



(3) 



(4) 



< J rnS m 



-S' + f7 m 6, t „&„ 



where a; and b m are bose operators and satisfy the 
relations [ai,a\] = 8u> and [b m , b^,] = S mm /. Here 
S' = S(ai) av = 5(1 — x) and x is the impurity con- 
centration. The bracket (• • -) a v means the average over 
all configurations. In order to introduce the decrease in 
the number of spin with impurity doping into our theory, 
we replace Sai with S(ai) av . When the spin 5 is sub- 
stituted by the nonmagnetic impurity, bonds with the 
exchange integral J are randomly configured. The ran- 
domness of both the on-site spin and the exchange inter- 
action must be considered. Since the Heisenberg model 
has only the interaction between spins, the randomness 
of the exchange integral can be treated in the frame of 
the coherent potential approximation (CPA) with the 
Tahir-Kheli's two-site approximation. And also the on- 
site randomness can be expressed by replacing 5 with 
5'. Rewriting the Hamiltonian (|l|) by eqs. (3) and (4) 
and ignoring terms more than the zeroth order of 1/S", 



we obtain 

H c = 25' J2 Jai°M a i + b fo + a ih + a \ h ]) ■ ( 5 ) 

<i,j> 

We introduce the Green's functions defined by 

Gf lr (t-t') = {{crmit); a v 4(t'W , (6) 
F^t-t , )^{{a m bl(ty,a v 4(t'W , (7) 

where ((A;B)) C = -i0(i - t')([A.B]) c and (A) c = 
Tr^4e _/3ff °/Tre^ /3ffC . Energy Fourier transforms of equa- 
tions of motion are expressed by 

luG c u ,(uj) = ai5 w + 5'y^ Jaiai+s 
s 

x + Fi+sA")) > ( 8 ) 

LOF^^Uj) = -5' ^ JCmVm+S 
S 

x (G c mW H + F^,H) , (9) 

where 5's are z vectors to nearest neighbors and G^, (uj) 
and (uj) are Fourier transforms of Gy/(t — t 1 ) and 
Fmi'(t ~ respectively. 

In order to treat the nonmagnetic impurity effect of the 
exchange integral, we introduce the coherent exchange 
integral J p (cl>).[15] Adding the same terms to both sides 
of eqs. (8) and (9), we can write eqs. (8) and (9) as a 
matrix expression: 



r-y = A c + V c g c 

- - i - f (r- 1 )/, (r- 1 ),, 

mj 

[uj - J p (uj)S'z}Su -J p (u>)S' J^S 
J p (uj)S'J2s s i,m+8 [uj + J p (u))S'z]6 jm 



(10) 



v c = 



v c V c ■ 

mi mj 



o 



,m+<5 

Using the relation r _1 r = 1, we obtain 



r = 



Tij 

nj 

(00)/, ,^ifc.(i-i) 



1 mi 1 mj 



(11) 



2 r (10)/ s ik-(m-i) 2 r (H) {, : \Jk(m-j) 



with 

rr o, H = -r 



(ii), 



-UJ 



uj + J p (uj)S'z 



.(01) 



^-[ Jp(w )5'zP(l- 7 2) 
-J p (ui)S'zj k 



, (12) 



^-[J p ( w )5'^(l- 7 2) 



. (13) 
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and 



Ik 



s 



(14) 



where k runs over a half of the original Brillouin zone. 
Introducing the T-matrix through the relation: 



y Cg c = T c TA c ^ 

from eqs. (10) and ( p^| ) we have 

g c = TA C + TTTA C 
T c = V C + V C TT C 



(15) 

(16) 
(17) 



We use the two-site approximation to solve eq. (17) for 
T c and obtain 



1 W 



01 



S'[JtJ (Tl - Jp(u)] 

i - S'[J(jq(ji - j p (^)]r 



where 



2 

N 



^ ,,,2 _ 



2J p (u)S'z(l - 7 2) 



ufi-[J p (u)S'z]*(l-rt) 



From eq. (18) the non-scattering condition (T c 
of CPA leads the coherent exchange integral as 



J p (w) 



2x 



i + r p (w) 



where 



N ijj 2 

k 



[j p (w)S'z]*(i->r*) 



(18) 



(19) 



= 



(20) 



(21) 



The average over impurity configurations is performed in 
the following way: the site is occupied by a magnetic 
atom and the conditional averaging is taken over the two 
occupational states of the site 1. [15] That is, the prob- 
ability that the site 1 is occupied by a magnetic atom, 
is 1 — x and the probability occupied by a nonmagnetic 
impurity is x. 

Averaging eq. (16) and assuming (T c A c ) av ~ 
(T c ) av (A c ) av , we obtain the Green's function: 

w(o-;)av + J p (uj)S'z 



(G? i M) av = - 



N ^ LU 2 

k 



[J p ( W )S'zP(l- 7 2) 



(22) 



If the Green's functions in eqs. (8) and (9) are expanded 
in a series of JaiUj, the corresponding term of the second 
term of the numerator of eq. ( p2[ ) is Ja^ajSz. This term 
does not change by multiplying Oi because af — Oj. In 
order to avoid the double count of the impurity effect, we 
replace J p (uj)S' z(ai) av with J p (uj)S i z. Thus we must pay 
particular attention to the average of products of 07 . For 
example, when V = I, (ai<Ji>) av = (ai) av (ai>) av = (<Ji)%, 
must be replaced by (cr/) av . The correlation function is 
obtained from the spectral relation: 



{crtaiai) = 



dcu- 



-Im<G^(o; + iO)) av . (23) 




Fig. 1. Impurity concentration dependence of the sublattice mag- 
netization m and the coherent exchange integral J p for S = 1/2. 
The critical impurity concentration x c is 0.303. The percolation 
threshold x p is 0.500. mo is the sublattice magnetization for a 
pure system. 



Since the spontaneous sublattice magnetization is ex- 
pressed by 



(<nS? 



= 5(1 -x)- (aiajai) , 



(24) 



where (A) = ((A) c ) av , the impurity concentration de- 
pendence of m is obtained from eqs. (|2C| ) - (24). 

We consider the case J p (ll> — 0). This case corresponds 
to the real excitation spectrum without the damping due 
to the scattering of spin-waves by the disorder. We ob- 
tain J p (0) from eqs. @ and ©: 

^p(0) 



= 1 - 2a; 



(25) 



at which J p is zero, is the 



J 

The concentration x p = 0.5 

percolation threshold. Substituting Jp(0) into eq. (|2 
we can analytically integrate eq. (g3|) and obtain the 
magnetization: 



m = {S+ 1 -){l-x) 



N ^ 
k 



1 



7 fc 



= -coth(^ J p S'z 
2 2 v 2 1 



ll) ■ (26) 



At T = we get 



m={S+\){\-x)-\ 



1 



dk 



(2tt)2 



(' 5 + 7T)( 1 - a; )- - 6966 2 



(27) 



The x dependence of J p and m are obtained from eqs. 
( p5| ) and (27) and are shown for 5 = 1/2 in Fig. 1. As x 
is increased for 5 = 1/2, both J p and m decrease and m 
vanishes at x c — 0.303. At this concentration, however, 
J p is still finite and J p vanishes at x p = 0.500. Thus 
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the state, with the infinite connecting percolation path 
and without LRO, exists in the concentration region of 
Xq X X ry . The critical concentrations of S — 1, 3/2, 
2, and 5/2 are 0.536, 0.652, 0.721 and 0.768, respectively, 
while the percolation thresholds of these spins are also 
0.500. Since x c > x p for S > 1, LRO vanishes at the per- 
colation threshold. The spin dependence of x c is caused 
by the decrease of the magnitude of the effective spin. 
On the other hand, x p is independent of S. And so, the 
disagreement between x c and x p appears. The obtained 
critical concentration x c — 0.303 for S = 1/2 is disagree- 
ment with x c = 0.07 obtained by the Kondo-Yamaji's 
method, [13] because the Kondo-Yamaji's method over- 
estimates the spin fluctuation. 

§3. Modified spin-wave theory of a 2D Heisen- 
berg antiferromagnet with nonmagnetic im- 
purities 

In order to investigate the state in the impurity con- 
centration region of x c < x < x p , we apply the modified 
spin- wave theory to the 2D antifcrromagnetic Heisenberg 
model with nonmagnetic impurities. We consider the fol- 
lowing Hamiltonian: 

H C = 2J ]T ^S^+M^a^-^T^S;), (28) 

<i,j> i&A j£B 

where [i is a Lagrange multiplier. The second term en- 
forces the constraint that the total staggered magnetiza- 
tion is zero. We introduce the antifcrromagnetic Dyson- 
Maleev transformation: 



get 



aiSf = S' - aiajai , aiS t 
criSf = (26" - aia\ai)aiai 

CfmSfn = — S' + <T m 6j ra 6 m , (J m S m 



(29) 



and obtain equations of motion of Green's functions as 

{u + n)G c w {u) =ai8 w (31) 

+ 2J ^ 0707+5 [S 1 ' — (<ri+8m+8bl +5 bi + s) c 
s 

+ (ai<Ji + saibi +s ) c }(G c u ,{uj) ~ Ff +S ^,{uj)) , 

(w-p^Cw) (32) 

= 2J cr m a m+ g[S' — (cr m +8 (J m+8 a l n+ s a m+s) c 
s 

+ (v m <J m +sal n+s bl n ) c }(G c m+S j,(u;) - F^ v {uS)) . 

We approximate terms in the bracket [• • •] in eqs. (31) 
and (32) with those averaged over impurity configura- 
tions, i.e., 

C(6) = S' - (<Ji + s<Ji + sb] +s bi+s) + (aiai +s aibi +s ) . (33) 

The correlation function C(S) must be obtained self- 
consistcntly. 

In much the same way as the above spin-wave theory, 
we introduce the coherent exchange integral </ P (w) and 



L g c = A c + V c g c 



(34) 



L = 



i = / (r- 1 ^ (r-% 



9 



(r )mi (r )mj 

[Lu + fi-2J p {Lu)J2 s C{d)]5 lt 
-2 J p (u) J2s C(5)S m+ Sa 

[u- Li + 2J p (aj)^ s C(5)]S mj 
Gw(u) \ A c = ( °ih> 

F mi'M ) ' \ 

v c V c ■ 

mi mj 

^Ejl 1 ' 17 ™^ — Jp(w)]C(5)8 m +S,i 

-2l28[J (T i cr i+s - Jp( w )]C(8)6i + 8,j 

-t^siJCrnOm+S ~ Jp(uj)]C(8)8 m j 

Using the coherent potential approximation with the 
two-site approximation, we obtain Green's functions and 
Jp(w), i.e., 

(G^M)av = IfX] 



where 



N 
2 



N — u 
k 



E: 



uj{(Jy) av + X n ik-(l-V) 
2 (1 - i 



■ > , ■> | > ■:> ' - (35) 



2-A2(l-ry2 7 2) 



Jk(m-l') 



(36) 



A = 2Jp(w)2;C(5) - ju , , 5 ^f (i) , (37) 



and 



with 



J 



2.r 



i + r P H ' 



t^ 2 - A 2 (l - T?) 
^c;2-A2(l-, ? 2 7 2) 



(38) 



(39) 



For simplicity, we consider the case that J p (u>) is real. 
For cases with the spin gap, J p (co) is real for all ener- 
gies within the gap. If the spin gap vanishes, then the 
dependence of J p {u>) on x is required to tend to that of 
the spin-wave theory. Therefore, we put lj = A^/l — rj = 
y/X(—fi). Then from eqs. (|38| ) and we obtain 



j P (VAra ) = l 
j 



(40) 



From the spectral relation (|23|) and eqs. (35), (36) and 
E^) correlation functions are obtained: 

(oia v a\a v ) = f(l - I') - -8 w {a v ) m , (41) 

(42) 



(<ricr m albl n ) = g(l - m) 
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where 



l coth( ^^l_^ 7 2 )e -ife-(/-I') 5 (43) 



2^,2 ^ e -ife- (J-m) _ ( 44 ) 



2 

X I CO th(^^l-^ 7 2) 

In the modified spin-wave theory we put the spontaneous 
sublattice magnetization zero. Since nonmagnetic im- 
purities with a completely random configuration cannot 
change rotational symmetry of the system, we can put 
the magnetization zero for impure systems: 

= faSF) = (S+i)fa) w -/(0). ( 45 ) 

From cqs. (41), (42) and ( f45| ) correlation functions of 
spin operators and C(8) are written by 

(a l a l ,S l -S l ,) = [f(l-l')} 2 for IV I, (46) 

(a l a m S l -S m ) = -[g(l-m)} 2 , (47) 

C(8) = g(S) . (48) 

The on-site correlation function at I' = I is (aiSi ■ Si) = 
S(S + l)(* l )„ = S(S + l)(l-x). 

At T = from eqs. (43) and (|||) we have 

1 



(49) 



We take the limit TV — > oo and rewrite the sum into the 
integral for the right-hand side of eq. ( |49| ) . Although the 
sum is divergent as r\ — > 1, the integral is not divergent. 
This contradiction can be avoided by introducing Bose- 
Einstein condensation. [16] Separating the divergent term 
at k = (0, 0) from the sum, we have 

1 f dk 1 

2.1 j^rJT 



(S+l)(i-x) = — i 



. (50) 

The mean-square root of the staggered magnetization m 
is defined by 



= ^£(-l) n (<W?o-Sn) 



(51) 



This is the long range order (LRO) parameter in the 
thermodynamic limit. From eqs. (43), (44), (46), (47) 
and (51) LRO parameter is given by 

2A 2 Y !-»7 2 7fc 2 



<oth 2 (§,/l-^ 7 2) 



4A 



(1-x). (52) 



For J p we take the limit T — > on finite lattice. 
From ?7 1 we have 

■,2-,2 



1 V l + r^7| 1 , 
2N 2 ^l-rf 1 l AN 

k K 



(53) 



In the limit N — > oo, integrating eq. ( p3|) except A; 
(0,0) and neglecting x/AN, we obtain 

1 1 



" N 2 \ ~ if (M) 

From eqs. ([30]) and ( |54] ) LRO parameter at T = is 
written by 

, I f dk 1 



7fc 



(5 + -)(!- a;)- 0.696602 



(55) 



Thus x c and greement with those of the spin- 

wave theory. Our modified spin-wave theory has the so- 
lution with rj 7^ 1 in x c < x < x p from eq. (fl9"|). The 
state in this concentration range is the disordered state 
with the spin gap. We conclude that La2Cui_ a; Mg :; .04 
has the disordered state with the spin gap in x c < x < 
x p , but K 2 Mni_ :E Mg\ I .F4 is in the Neel state at all con- 
centrations x < x p . This disordered state is presumably 
due to the construction of the singlet dimer bonds. 

§4. Summary and Discussion 

In the present paper we have applied the spin-wave 
theory and the coherent potential approximation (CPA) 
to the impure 2D Heisenberg antiferromagnet with spin 
S. As magnetic ions of spin S are substituted by non- 
magnetic impurities, both the number of magnetic ions 
and interacting spin pairs are decreased. We take these 
two effects into account by substituting 5(1— x) for S and 
using CPA to the exchange integral JoiOj. At T = the 
critical concentration x c is 0.303 for 5=1/2 and 0.500 
for S > 1. On the other hand, the percolation thresh- 
old x p is 0.500 for all spin S's. Thus for S = 1/2 the 
state, with the infinite connecting percolation path and 
without long range order, exists in x c < x < x p . For 
S > 1 the long range order vanishes at the percolation 
threshold. The decrease in the effective spin leads the 
difference between x c and x p for S = 1/2. To investigate 
the state in x c < x < x p , we have studied by the modified 
spin- wave theory. As a result, this state is the disordered 
state characterized by the finite spin gap and the expo- 
nentially decaying correlation function. This disordered 
state is presumably due to the construction of the singlet 
dimer bonds. 

Thus our results qualitatively explain the experimen- 
tal results: the critical impurity concentration x c for 
the Neel state of La2Cui_ a: Mg !( .04 (S — 1/2) is smaller 
than the percolation threshold x p = 0.41, on the other 
hand, x c of KaMni^Mg^ (S = 5/2) is equal to 
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x p . La 2 Cui_ :E Mg a ,04 would be the disordered state in 
x c < x < x p . We expect to experimentally observe the 
spin gap in low temperatures for La2Cui_ a; Mg a ,04 in 
x c < x < Xp or for other materials with S = 1/2. 
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